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Abstract 



We address the fi problem of gauge mediation by considering a singlet chiral superfield 
coupled to the Higgs and messenger fields. We compute the soft terms generated below 
the messenger scale and study the phenomenological consequences of the model. The 
experimental bound on the Higgs mass provides a severe constraint that identifies three 
special regions of parameters where the mass spectrum and the collider signatures can 
be distinct from ordinary gauge mediation. 



1 The Origin of /i in Gauge Mediation 



Gauge mediation [H [21 [3l H] is one of the most appealing realizations of the supersymmetric 
extensions of the Standard Model (SM), because of its high predictive power and ultraviolet 
insensitivity. Indeed, the induced soft terms are computable in terms of few parameters (gen- 
erally three) and do not lead to unacceptably large flavor violations, under mild assumptions 
on the unknown interactions responsible for generating Yukawa couplings. 

There are two ingredients of gauge mediation that are still obscure. One is the seed of 
supersymmetry breaking, which is expected to have a dynamical origin, in order to explain 
naturally the emergence of mass scales much smaller than the Planck mass Mpi. New in- 
teractions have to transfer the original supersymmetry breaking from the hidden sector to 
the messenger fields. Attempts to simplify this structure and unify the hidden and mes- 
senger sectors have faced various difficulties. One of the problems is that, once a hidden 
sector with supersymmetry breaking is found, the couplings to the messengers allow for 
new supersymmetric vacua. To address this problem, models were constructed E] around 
metastable false vacua where supersymmetry is broken, although the global minimum re- 
mains supersymmetric. More recently, a new framework of theories with these properties 
has been discovered [7|, showing that this approach is fairly general. This new development 
has revived the interest in gauge mediation and led to the construction of many interesting 
models [8]. 

The second obscure ingredient is the seed of Peccei-Quinn (PQ) symmetry breaking or, in 
other words, the origin of /i, the higgsino mass, and 5^, the square mass mixing the two scalar 
Higgs doublets Hd and Hu- The first aspect of this problem, common to all supersymmetric 
models, is how to relate /i to the soft masses of the other supersymmetric particles. This 
problem is solved by assuming that the PQ symmetry is exact in the supersymmetric limit, 
while /i is induced by supersymmetry-breaking effects [9j. The second aspect appears only in 
theories with computable soft terms, like gauge mediation, and it is expressed by the generic 
prediction 

where F represents the value of the hidden-sector auxiliary field, and M is the mass of the 
mediating field. Equation ([1]) is the consequence of a generic coupling in the Kahler potential 
between the Higgs bilinear HdHu and a spurionic superfield X = 1 + 9'^F/M 



where a represents the product of coupling constants and possible loop factors required to 




(2) 
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generate the effective interaction. The operator in eq. ([2]) simultaneously generates both /i 
and Bfj^, leading to eq. independently of the actual value of the loop-suppression factor 
a. In gravity-mediated models, F/M corresponds to the natural scale of soft terms, and 
eq. ([1]) is fully satisfactory. However, in theories where the soft terms are derived from F/M 
through computable loop effects, like gauge or anomaly [10] or gaugino [TT] mediation, eq. ([T]) 
predicts that the ratio B^/ fi is parametrically too large, requiring an unnatural fine tuning. 

One solution [12] is to construct models where operators of the form ([2]) are absent, while 
couplings to the hidden sector generate only the structure 

J d'eH,H^D'f{X,X^), (3) 

where is the supersymmetric covariant derivative. Since is an antichiral 

superfield, the operator in eq. ([3]) generates /i but not B^, which is then induced at a 
higher order in perturbation theory. Other solutions use other dynamical scales present in 
the hidden sector [13] or required by the cancellation of the cosmological constant [H] to 
reproduce acceptable values of /i and B^. It is also possible to construct models with flavor 
symmetries [15] or an R-symmetry [16], leading to selection rules that invalidate eq. ([1]). 

An alternative approach is to introduce in the low-energy theory a new SM singlet field 
coupled to the Higgs bilinear in the superpotential 

W = \NH,H^ - ^N^ (4) 

A symmetry forbids a bare fi term, and the coupling k is needed to break the global 
PQ symmetry. The effective /i and 5^ terms can now be entirely generated by low-energy 
dynamics and /i = X{N), B^ = A(F/v) ~ (A^)^, circumventing eq. ([1]). 

In theories like gravity mediation, where /i is correctly generated by supersymmetry 
breaking and eq. ([T]) is successful, the introduction of the singlet A^ does not appear to be 
well motivated. Not only is it superfluous, but it also introduces a proliferation of new 
unknown parameters in the soft terms. Moreover, a light singlet can potentially destabilize 
the hierarchy [T7], as is the case when we embed the superpotential in eq. (j4]) into a GUT. 
The situation is quite different in a theory with low supersymmetry-breaking scale and 
computable soft terms, like gauge mediation. In this case, the introduction of A^ is essential 
to bypass eq. ([1]) and, in principle, it can be done at the price of only two new parameters 
A and k (to be compared with /i and -B^ of the minimal supersymmetric SM) in the full 
Lagrangian, including soft terms. Moreover, as long as y/F < 10^ GeV [4], the coupling 
of A^ to GUT fields does not destabilize the hierarchy [18], and it can even be used for 
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the sliding-singlet mechanism jlH] to explain the Higgs doublet-triplet splitting, in the limit 
A; -> pn]. 

The use of to generate the fi term in gauge mediation was immediately suggested in 
the original paper on the subject pLj, but it was also found that the specific form of the soft 
terms in gauge mediation does not allow for a correct pattern of electroweak breaking with 
an acceptable mass spectrum. The main difficulty lies in generating a sufficiently large value 
of {N), which requires either a negative soft square mass for the scalar field A^, or large 
A-terms for A and k interactions. Neither of these conditions can be satisfactorily obtained 
in gauge mediation, unless one introduces new hght fields coupled to P, multiple singlets 
with appropriately adjusted couplings [2], higher-dimensional interactions of with specific 
values of the exponents [3], or modifies the theory to include a new U(l) gauge group, under 
which A^ is charged, with new associated fields [21] . A thorough analysis of these possibilities 
has been presented in ref. [22] . 

In ref. |2^ it was pointed out that a negative square mass for A^ and non-vanishing 
trilinears can be obtained if the singlet is directly coupled to the messenger fields $ in the 
superpotential 

W = X {ki^I^I + K2^2'^2) +^N^l^2, (5) 

where X is the hidden-sector superfield containing the Goldstino. The form of eq. ([5]) 
can be guaranteed by symmetries, e.g. by a discrete with Z3[$i] = 2'3[<l>2] = —1/3, 
■^3[^2] = ^ai^*!] = ^3[A^] = 1/3, Z^lX] = 0, broken only at the weak scale. The doubling of 
the messenger field is necessary to avoid a kinetic mixing between X and A^. Indeed, if both 
X and A^ coupled to the same bilinear $$, below the messenger mass M we would find the 
one-loop mixing in the effective Kahler potential 

where (i$ is the dimensionality of the gauge representation of $. This generates a tadpole 
for the scalar field A^ {Vcs = (^(i$/167r^)A^F^/M) that destabilizes the weak scale, unless 
a/F < TeV / v^- Therefore the doubling of messengers is necessary, unless we accept the 
unnatural case of an extremely small value of ^ [23] . Instead of coupling the singlet to both 
messengers as in eq. ([5]), one could also consider coupling it to (the weak-doublet component 
of) just one messenger, as in N^iHu or NHd^2- This alternative was discussed in ref. [25] . 

In this paper we study the structure and the phenomenology of gauge mediation with an 
extra singlet A^ coupled to Higgs and messenger fields as in eqs. (jlj) and (j5]). This variation 
of the ordinary SM extension with gauge-mediated supersymmetry breaking (GMSB) will 



3 



be referred to as the N-GMSB model. In section 2 we compute the soft terms induced in 
the effective theory below the messenger scale. In sections 3 and 4 we study the vacuum 
structure and the phenomenology of N-GMSB. Our results are summarized in section 5. 



2 Generating the soft terms 

We now want to compute the soft terms of N-GMSB, the supersymmetric SM with gauge 
mediation augmented by the superpotential interactions in eqs. (jll) and ([5]). We start by 
treating X = M + as a background non-dynamical field, and we will later comment on 
the case in which X can propagate. We also assume that $ ($) belongs to a fundamental 
(antifundamental) of SU{b) and we introduce separate couplings for the interactions of 
the SU{2) doublet ($^) and the SU{?>) triplet (<l>^) contained in the messenger multiplet, 
expanding the relevant terms in the superpotential as 

A k 
W = XY, + f^^^^^J) + N (eD$f + ir^l^l) + XNH.Hu - -N\ (7) 

i=i 

We define the supersymmetry-breaking mass and interaction terms for the Higgs bosons 
of the effective theory valid below the messenger mass as 

Koft = rhl^\Hu? + m^Ji^dP + m^|iV|' + [xA^NH^H^ - ^A^iV^ + h.c.^ . (8) 

Soft terms are obtained by integrating out the messengers at one loop, for A terms, and at 
two loops, for scalar masses. Instead of computing the full set of diagrams, a rather daunting 
task, we use the method proposed in refs. [231 EE] to extract supersymmetry-breaking effects 
from wave-function renormalization. 

We first write the one-loop Renormalization Group (RG) equations for the field wave- 
function renormalizations and the couphng constants Aj as 

d\n.Zn 



dlnQ 



7„ a = N, Hd, (9) 



^j^g =/5a,, \ = \,k,iD,iT,g,g',gs- (10) 

Here Q is the renormalization scale, 7^ are the anomalous dimensions, and I3\. are the beta 
functions. Retaining only the leading terms of an expansion in powers of F/M, the soft 
supersymmetry-breaking parameters in eq. ([8]) are given by (see appendix B) 
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where we have defined AX = [X'^^^ —X^~^]q=m (with X = /5a,, 7a) as the discontinuity at the 
messenger scale, and X^^^ are the values of X in the theory above and below M, respectively. 
Z'^ is proportional to the discontinuity of the anomalous dimension at the messenger scale 
M, and Z'^ depends on a combination of the discontinuities of the anomalous dimension and 
beta-functions. Such discontinuities can be present if some particles contributing to Px. or 
7 are integrated out at the scale M. 

To obtain explicit formulae for the matching conditions on the soft terms at the scale M 
we only need to specify the anomalous dimensions, beta functions and their discontinuities 
at M. This can be easily done from the RG equations given in appendix A. The anomalous 
dimensions and their discontinuities are given by 
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»7r^ ~ ' 87r" 

In eq. ( IT4l) . we have not specified the "Yukawa" contribution, which is different for Hd 
and Hu, since it does not lead to any discontinuity at the scale M and therefore does not 
contribute to soft masses up to two-loop order. The beta-functions and their non-vanishing 
discontinuities that contribute to Z'^ in eq. ( JT3ll are 
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Here the gauge couplings gi are ordered as {g',g,gs) and the constants q are (5/3, 1, 1); n 
is the number of messenger pairs (we take n = 2). Finally, the matching conditions on the 
soft terms at the scale M are explicitly written as 
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Here all couplings and parameters are evaluated at Q = M. 

Note that there is no one-loop contribution to of order F'^/M'^. Indeed, the messenger 
interactions in eq. ([7]) are invariant under independent chiral reparametrizations of the fields 
N and X (with messenger fields transforming appropriately), therefore the one- loop Kahler 
potential must be of the form J d'^ON'^N In X'^X and cannot induce a soft mass for A^. 

A one-loop contribution to can be generated only at higher orders in the F/M'^ 
expansion, and we find 



m = n Z ■ (21) 



1 — + 2x\nx 

This contribution is always negative. However, it is negligible with respect to the one in 
eq. ( |T9l) . as long as M ^ AirF/M ^ 10^ TeV. As we will see in section HJ this condition 
is satisfied in our study, since the Higgs mass bound selects large values of M. A variant 
to induce a one- loop contribution to m% of order F^/M"^ is given by models with several 
hidden-sector fields with non- vanishing vacuum expectation values (vevs). This case can be 
parametrized by a superpotential interaction + X2^2^2, with Xi 2 = Mi 2 + ^^-^1.2 

where M12 and Fi 2 are independent. Now we have to consider two different messenger 
thresholds and the one-loop contribution is given by 



w \m, mJ )' 



= (1 - ^) + (1 + ^) 1^^] • (23) 



Note that the contribution in eq. (]22i) is always positive, and vanishes when the supersym- 
metry breaking is universal, Fi/Mi = F2/M2, (as is the case for a single X field, considered 
in this paper, where the coupling constants /ti,2 drop out from the ratio -Fi, 2/^1^2) or when 
one messenger threshold decouples (Mi or M2 00). 

The matching conditions at the messenger scale on the soft masses for gauginos, squarks 
and sleptons are given by the usual expressions of gauge mediation 

M, = nc,^^, (24) 
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where the coefficients Cj and n are given below eq. f|T6|) . and C/ is the quadratic Casimir 
invariant for the scalar / under the gauge group with coupling ctj. The matching conditions 
at the messenger scale on the trilinear A-terms corresponding to Yukawa interactions vanish 
at leading order, while nonzero values are generated at the weak scale by RG evolution. 

As is well known, the couplings ki^2 of X to the messenger fields do not affect the soft 
terms, since they drop out of the ratio F/M. However, when messengers are coupled to 
N, a propagating X field gives a two-loop diagram that contributes to m^. Indeed, with a 
dynamical X we obtain an extra contribution to /^^^''^ giving 



+ K 



This leads to an extra term to be added to eq. f|T9|) 

1 
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In the rest of the paper we will restrict our analysis to the case in which X is a spurion 
representing only the mass parameters M and F (or to the case in which X propagates, but 



are negligible with respect to the other coupling constants) and neglect the contribution 
in eq. 



3 Vacuum Structure and Higgs Boson Masses 

To determine the mass spectrum of the low-energy limit of N-GMSB, which essentially can 
be viewed as a constrained version of the Next-to- Minimal Supersymmetric Standard Model 
(N-MSSM), we must compute all the Lagrangian parameters at some renormalization scale 
of the order of the weak scale, where we impose the minimization conditions of the Higgs 
potential. The model has five unknown input parameters: the singlet couplings A and fc; 
the messenger mass M; the effective supersymmetry breaking scale F jM\ the unified value 
^jj = ^r,z)(MGUT) for the singlet-messenger couplings at the GUT scale (defined as the 
scale where the couplings g and ^^/5/3g' meet). Other required inputs are the gauge and 
third-family Yukawa couplings, which we extract at a low reference scale equal to the pole 
top mass Mt = 170.9 GeV [27] from the known values [28] of the fermion masses and of 
the SM input parameters Gp, Mz, sin^ and as(M^) . We use tree-level formulae for the 
determination of all the couplings but the top Yukawa coupling ht, for which we include 
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one-loop corrections. The soft supersymmetry-breaking masses and interaction terms for 
Higgs bosons, gauginos and sfermions are determined at the messenger scale M by means of 
eqs. f|T7j) - f|T9l) and eqs. fl2^ - fl25|) . Finally, we determine all the parameters of the N-GMSB 
Lagrangian at a renormalization scale Ms that, in order to minimize the dominant 0{hf) 
one-loop corrections to the Higgs potential, we choose as the geometric average of the two 
stop masses, i.e. Ms = ^JmijrTi^. To this purpose we use the RG equations of the effective 
theories valid between the different mass scales: SM between Mt and Mg; N-MSSM between 
Ms and M; N-GMSB including the messenger sectoJl] between M and Mqut- The explicit 
formulae for the RG equations are given in the appendix A. Since the boundary conditions 
on the various parameters are given at different renormalization scales, and some of them 
depend on the vevs {Hu) and (i/^) determined by the minimization of the Higgs potential, 
we need to iterate the procedure until it converges. 

The tree-level scalar potential along the neutral components of the fields and N is 



\Hd\ 



+ 



k 



XA^NH.H^ - -AkN' + h.c. ) + mj^jH^l'' + mj^jH,]'' + m%\N\\ {2i 



The minimization conditions of the scalar potential with respect to the three Higgs fields 
allow us to determine the vevs , {Ha) and (A^). In practice, we treat the electroweak 
symmetry-breaking scale = (Hu)"^ + (Hj)'^ ^ (174 GeV)^ as an input parameter extracted 
at Q = Mt from the Fermi constant Gp and evolved up to Q = Ms with the SM RG 
equations. The minimization conditions can therefore be used to determine one of the 
unknown input parameters, reducing their number to four. In terms of parameters computed 
at the scale Ms, the minimization conditions can be expressed as 



^h,-^hJ^^ P g^ + g 
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tan2/5 
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k\ sin 2/3 X'^v^ sin"^ 2p' 



(30) 



(31) 



where, to highlight the analogy between eqs. (I29l)-(l30l) and the corresponding minimization 
conditions in the usual MSSM, we define tan/3 = (Hy) / {H^) and introduce the quantities /x 



^We neglect possible self-interactions in the hidden sector. Otherwise, as argued in ref. 
consider one more effective theory valid between the scales M and 



we should 
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and Bfj^ : 

fi = X{N), = jfi''-Axfi-^sm2P . (32) 

Eqs. (l29l) - fl3Tl) depend non-trivially on the various parameters and must be solved numeri- 
cally. They provide us with values a.t Q = Ms for (A^), tan/5 and a third parameter that 
we choose to be k. The remaining input parameters are thus M, F/M, and A (the latter 
given at the scale Ms). Without loss of generality, we can take A real and positive and exploit 
the freedom to redefine the phases of the fields to make sure that tan f3 is also positive; this 
involves flipping the signs of (A^) and k if the numerical solution of the minimum equations 
gives a negative value for tan/5. We also choose a basis in which F/M is real and positive, 
so that the gaugino masses in eq. are positive. In the general N-MSSM, the relative 
phases between the A-terms and the gaugino masses cannot be removed and are physical 
sources of CP violation. However, in the N-GMSB these phases are all zero and, in the field 
basis we have chosen. Ax and Ak at the messenger scale turn out to be real and negative, 
see eq. ( IT71) . 

Before discussing our treatment of the radiative corrections in the N-MSSM Higgs sector 
and moving on to the numerical analysis, we present some analytical considerations that 
help understanding the vacuum structure of the theory. After LEP unsuccessful searches for 
the Higgs boson and for new particles, supersymmetric models suffer from a mild fine-tuning 
problem that requires a certain separation of scales between v and the superparticle masses. 
Therefore the only acceptable region of parameters has to lie very close to the "critical line" 
separating the phases with broken and unbroken electroweak symmetry [30] . In practice, 
this means that we can find (A^) by setting f = in eq. fl3Tl) . and then imposing the critical 
condition for electroweak breaking on the effective Higgs potential at fixed A^ background 
value. From eq. (j3Tl) we obtain 

^1 = -AkW + 0{v^) , w = , z=-^. (33) 

k 

The non-trivial vacuum for A^, corresponding to eq. (1331) . exists only for z < 1/8, but we 
have to impose z < 1/9 to insure that this vacuum is deeper than the origin (A^) = 0. This 
condition then leads to w > 1/3. 

Equations (129|) and (130|) . in the limit v 0, correspond to the critical condition that 
the origin of the effective Higgs potential with (A^) fixed has locally one flat direction and 
non-negative second derivatives: 

+ /i') + = Bl , (34) 

^L + ^L+V>0. (35) 
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Combining eqs. 
\-k plane: 



and neglecting terms of 0{v'^) we obtain the critical line in the 
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A^w^- ^^^^ 

Equation (136|1 shows that the critical line gives an approximately linear relation between A 
and k, distorted only by the small mfj^ contribution and by RG effects. At large values of A 
(and k), the critical line is interrupted either by perturbative constraints on A and the top 
Yukawa coupling, or by the appearance of a minimum with {Hu) 7^ and (Hd) = {N) = 0, 
which becomes deeper than the correct vacuum unless 



■ w 



1 

"'-3 



(37) 



It is also interesting to note that the combination of the minimization conditions of the 
scalar potential with the boundary conditions on the Higgs trilinear couplings leads to a 
definite prediction for the sign[^ of /i. Indeed, it can be seen from eqs. (|30|) and fl35l) that the 
condition tan/3 > requires to be positive. Combining eqs. (132|) and (l33ll . and neglecting 
terms of 0{v^), one gets 

B^c:^fi{AkW-Ax) . (38) 



At the messenger scale Ak is negative and equal to 3^4^, see eq. ( 1171) . If the effect of 
the RG evolution of the soft supersymmetry-breaking parameters down to the scale Ms 
is neglected, the condition w > 1/3 constrains ^ to be always negative. In practice we 
find that, even though the RG evolution can alter the relation between Ak and A\, all the 
phenomenologically viable solutions to the minimization conditions of the scalar potential 
have indeed fi < 0. 

In the limit (A^) ^ v the tree-level squared masses of the two CP-odd and three CP-even 
neutral scalars are 
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(39) 



(40) 



(41) 



^ Our choice for the sign of the superpotential term XNHdHu in eq. ([7]) corresponds to the convention in 
which the off-diagonal element of the stop mass matrix contains mt /i cot (3 and the chargino mass matrix 
contains — /i . 
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Analogously to the decoupling limit of the usual MSSM, the lightest CP-even Higgs boson 
hi has SM-like couplings to fermions and gauge bosons, and its mass is of O(f^), with an 
additional contribution - in curly brackets in eq. fj40l) - that is not present in the MSSM; the 
condition w > 1/3 also ensures that the second term in the curly brackets is always negative. 
The CP-even boson /12 and the CP-odd boson ai are heavy and have couplings similar to 
those of the MSSM Higgs bosons H and A (the same applies to the charged boson). The 
CP-even boson /13 and the CP-odd boson 02 are mostly singlet and are mostly decoupled 
from matter fields. We find > for z < —1, and m/13 < for —l<z< 1/9. 

Detailed studies of the N-MSSM Higgs sector date back to the nineties [31]. It is also 
well known that in supersymmetric models the radiative corrections involving top and stop 
loops can give a substantial contribution to the Higgs boson masses [32], and they must be 
taken into account for a meaningful comparison with the mass bounds from direct searches 
at LEP [33] . The dominant one- loop corrections, enhanced by four powers of the top Yukawa 
coupling ht, can be computed in the effective potential approach. The radiatively corrected 
effective potential for the Higgs fields can be written as Ves = Vq + AV, where Vq is given in 
eq. (l28ll and the correction AV is expressed in terms of field- dependent masses and mixing 
angles. The radiative corrections to the minimization conditions of the scalar potential are 
taken into account by replacing in eqs. (l29|) - (l3T]l 

.2^1 dAV 



, . , (H,, H^, N) , (42) 

where the subscript "min" means that the Higgs fields are set to their vev after computing 
the derivative of the potential. The radiative corrections to the 3x3 mass matrices for the 
CP-even and CP-odd Higgs bosons are in turn: 

1 d'^AV 



^ ^''y ~ 2 9Re4<9Re< 



mm 



2 dim (hi d Im 1 



(43) 



We have explicitly computed the dominant 0{hf) corrections to the minimization condi- 
tions and to the Higgs mass matrices given in eqs. (H2l) and fH3l) . and checked that our results 
agree with those available in the literature [M]. In addition, we include in our determina- 
tion of the Higgs masses the one-loop leading logarithmic corrections of 0{h^g'^, h'^\^) (by 
multiplying the mass matrices by appropriate wave-function-renormalization factors) and 
the two-loop leading logarithmic corrections of 0{h^g1, h^). Finally, after diagonalizing the 
Higgs mass matrices we include the one-loop leading logarithmic corrections of (9(A^) to the 
mass of the lightest CP-even Higgs boson hi, computed in the limit where (A^) ^ v. These 
corrections are accounted for by the term 

A<--^ln^, (44) 
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where the hghtest-Higgs quartic couphng is defined as Xh-^ = m|y(2f^), and the tree- level 
mass of hi in the limit (A^) ^ v was given in eq. (HOl) . The correction Am|^ in eq. (jSj) can 
be numerically relevant only if A is fairly large. It also includes some (not all) of the one-loop 
leading logarithmic corrections that involve the electroweak couplings, but the contribution 
of such terms is generally small. 

We have compared the results of our calculation of the Higgs boson masses with those of 
the general N-MSSM model using the public computer code NMHDECAY [35], which includes 
also the corrections controlled by the bottom Yukawa coupling as well as a more refined 
treatment of the one-loop leading logarithmic corrections controlled by powers of A and 
of the electroweak couplings. We find that, in points of the parameter space that will be 
relevant to our analysis, the two determinations of the lightest CP-even Higgs mass m/^^ 
agree within 5 GeV, with NMHDECAY predicting in general smaller values of m/j^ than our 
calculation. We consider this agreement satisfactory, given the approximations involved in 
our calculation - we neglect the one-loop electroweak corrections and the two-loop non- 
leading-logarithmic terms of 0{hfgg, hf) - and the unavoidable uncertainty coming from 
uncomputed higher-order corrections. 

4 Phenomenology 

As discussed in the previous section, the requirement of successful breaking of the electroweak 
symmetry reduces the free parameters of the model to four: the messenger mass M, the 
effective supersymmetry-breaking scale -F/M, the GUT-scale singlet-messenger coupling E,u 
and the singlet-Higgs coupling A computed at a renormalization scale Ms of the order of the 
average stop mass. The parameters M and F/M should be chosen in such a way that the 
0{h1) radiative corrections involving top and stop loops are large enough to lift m/u above 
the bound from direct searches at LEP (large values of the supersymmetric scale also imply 
that the heavy Higgs bosons are essentially decoupled, thus the LEP lower bound of 114.4 
GeV |33] on the mass of a SM-like Higgs boson applies). In the usual GMSB the condition 
that the trilinear Higgs-stop coupling At be zero at the messenger scale results in a small 
stop mixing at the weak scale. Therefore, a large value of M5, greater than a (few) TeV, is 
required to make large enough. In the model with an additional singlet, on the other 
hand, positive contributions to can arise when A is large and tan/3 is small, see eq. ( HOj) . 
However, the conditions of correct electroweak symmetry breaking and perturbativity of the 
couplings up to the GUT scale require \{Ms) ^ 0.55, and a sizeable contribution to m/^^ 
from radiative corrections remains necessary. We will therefore choose a value of F/M large 
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enough to result in an average stop mass of the order of 2 TeV. A large messenger mass M 
is also required to ensure that a sizeable value of At is generated by the RG evolution down 
to the weak scale. 



Figure [T] shows the values of tan (3 that result from the minimization of the scalar potential 
over the plane - K^s)- We choose M = 10^^ GeV and F/M = 1.72 x 10^ GeV, resuhing 



allowed region, which is delimited by the solid (red) line, no satisfactory solution to the 
minimization conditions of the scalar potential is found. In particular, the points on the left 
of (and above) the allowed region are ruled out because the minimum with (Hu) ^ and 
{Hfi) = (N) = is deeper than the correct vacuum, see eq. (l37j) . The points on the right of 
the allowed region for X{Ms) ^ 0.5 are ruled out by the requirement that the top Yukawa 
coupling be perturbative up to the GUT scale. For \{Ms) ^ 0.5 we find tan/3 > 1.5, 
resulting in a not-too-large top Yukawa coupling. Therefore, the upper-right arm of the 
allowed region can extend up to ~ 1.1, where the couplings A and k approach the 
perturbativity bound at the GUT scale. 

For values of C,u in the vicinity of the right edge of the allowed region tan (3 goes down to 
about 1.4. When C,u decreases, tan/5 increases, reaching values greater than 10 in the leftmost 
corner of the allowed region. In addition, values of tan/5 greater than 10 are obtained for 
large C,u and small A. This behaviour can be qualitatively understood by considering that 
eqs. (l30l) and (I32 l) -( l33l) . in the limit of large tan/?, reduce to 



For small A and small ^[/, k/X aX the weak scale is determined by the critical line in eq. (1361) 
to be small (and w ^ 1/3). On the other hand, for large C,u we approach the condition 
that the vacuum with non- vanishing (A^) is nearly degenerate with the origin, i.e. w ~ 1/3. 
Since the boundary condition at the messenger scale in eq. f|T71) gives Ak = 3A\, there is an 
approximate cancellation in eq. (1451) . 

The values of /i that result from the minimization of the scalar potential are in general of 
the order of the stop masses, but they are inversely correlated to the values of tan jS shown 
in fig. [H In the region where tan/? is small /i gets as large as 2.8 TeV, while in the regions 
where tan/? is large /i goes down to 1.4 TeV. Indeed, one can see from eq. (|29l) that a value 
of tan j3 close to 1 enhances /i, both because of the factor tan^ /5 — 1 in the denominator and 
because a smaller tan/? results in a larger ht, enhancing the stop contribution to the running 



in an average stop mass Ms = ^ym^^m^ that varies between 1.9 and 2.1 TeV. Outside the 




(45) 
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Figure 2: Mass of the lightest CP-even Higgs boson hi in the ^u~^{Ms) plane, for M = 10 
GeV and F/M = 1.72 x 10^ GeV. 
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Figure [2] shows the mass of the hghtest CP-even Higgs boson hi over the plane ^u — X{Ms), 
with the same inputs as in fig. [H It can be seen from fig. [2] that there are three separate 
regions in the plane — X{Ms) where the mass of the lightest CP-even Higgs boson is 
sufficiently large. These three regions can be characterized as follows: 

Region I. In the lower-left corner of the plot, where both A and C,u ^ the singlet vev 
(N) is generated by a large and negative value of m%, and tan /3 is proportional to X/k, taking 
on relatively large values 4 ^ tan/? ^ 20. The tree-level mass of the lightest Higgs boson 
is dominated by the first term in the r.h.s. of eq. (HQ]), as in the MSSM, and the radiative 
corrections lift m/j^ above the LEP bound. However, for extremely small values of the 
ratio X/k becomes very large, and the term unsuppressed by sin 2/? in the curly brackets of 
eq. ( I40l) gives a large and negative contribution to m/^, dragging it again below the LEP 
bound. This is one of the reasons why we cannot consider the N-GMSB with C,u = 0. The 
masses of the MSSM-like heavy Higgs bosons /i2 and ai are of the order of /i, while and 02 
become much lighter as tan/3 grows, since m/13 ~ 2 (k/X) jj, and rriaj ~ \/sJw (k/X) jj,. Note 
that 02 becomes an approximate i?-axion, because Ak 0. The fermionic component of 
(singlino) has mass ~ '^hg- Since in this region k/X is small, the singlet-like scalars and 
the singlino are considerably lighter than the other non-SM particles, and the singlino can 
be the NLSP. The light pseudoscalar has sufficiently small couplings to escape LEP bounds. 
The NNLSP is a bino-like neutralino, which decays into the singlino and a (real or virtual) 
SM-like Higgs boson hi. In this scenario, the decay chains of supersymmetric particles end 
with the NLSP singlino decaying into a pseudoscalar singlet and a gravitino, with a rate 



This process is a peculiar characteristic of the N-GMSB model. However, the NLSP decay 



by the LEP bound on m/^^. 

Region II. In the lower-right corner of the plot, where A is small but is large, the soft 
mass is positive, and the vev (A^) is generated by the large value of Al, see eqs. ( fTTl) and 
(fT9|) . The parameter tan/5 is large due to the cancellation in eq. (H5|) obtained for w ~ 1/3, 
and, in contrast with what happens in region I, the ratio A;/A is large. As a result, the 
negative and tan /3-unsuppressed contributions to in eq. (HOi) are not important, and the 
tree-level mass of the lightest Higgs boson is approximately equal to Mz (as in the MSSM). 





(46) 



16 vr Mi F2 



can occur inside the detector only for ^/F roughly smaller than 10^ GeV, a region disfavored 
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Figure 3: Upper bound on the mass of the hghtest CP-even Higgs boson hi in N-GMSB and 
GMSB as a function of tan/?, for M = lO^^ GeV and F/M = 1.72 x 10^ GeV. 

The other particles have masses maj/S ~ ~ m/13 ~ {k/X)fi, while h2 and Oi have 

masses of order fi. Due to the large value of k/X the singlet-like scalars and the singlino are 
much heavier than the other scalars and neutralinos, making the particle spectrum similar 
to the one of ordinary gauge mediation. 

Region III. The last region with relatively large lies at large values of A and close to 
the right edge of the region allowed by perturbativity of the couplings, where tan (3 < 2. The 
soft parameters and are large and negative, and they both contribute to generating 
(A^). The ratio k/X is close to 1, therefore all the heavy scalars, as well as the higgsinos 
and the singlino, have masses of the order of fi, while the NLSP is the bino-like neutralino. 
Concerning the mass of the lightest Higgs boson hi, the first term in the r.h.s. of eq. (BUI) 
is suppressed by the low value of tan /5, but the term f ^ sin^ 2/3 is sizeable and hfts m/^ 
above the LEP bound. 

If we give up the requirement that the couplings be perturbative up to the GUT scale, 
considering only their evolution up to a relatively small messenger scale, we can accommodate 
larger values of X{Ms), resulting in a larger tree-level contribution to rrih^. For example, for 
M = 10^ GeV, F/M = 1.5 x 10^ GeV (so that the stop masses are of the order of 2 TeV) 
and X{Ms) ~ 0.7 we can find a range of values of ^i:),t{Ms) for which ~ 150 GeV. 
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It is interesting to compare the results for obtained in the N-GMSB with those that 
are obtained in the usual GMSB for the same values of M, F/M, and tan/5. The result is 
summarized in fig. [3l where we show the maximal value of rrih-^ as a function of tan (3 for the 
GMSB and for the three phenomenologically viable regions of the N-GMSB. In region III, 
for low tan (3, the Higgs mass can be larger than the corresponding value in GMSB. Region 
II gives a prediction for the maximum value of m/^^ that is identical to the one of GMSB. 
Finally, the upper bound on m/j^ in region I is approximately equal to the one of GMSB 
for intermediate values of tan j3, but becomes smaller at large tan (3, because of the negative 
contribution —X^v'^/k'^ in eq. ( HOl) . 

5 Conclusions 

There are two aspects of the /i problem. The first is related to the absence of a /i term in 
the limit of exact supersymmetry and to its generation from supersymmetry breaking. The 
second aspect (the "S^ problem") is related to the generic expectation B^/n ~ F/M, and 
it is present only in models where the soft terms are calculable and turn out parametrically 
smaller than the original seed of supersymmetry breaking F/M. Gauge mediation belongs to 
this class of models. For these models, as opposed to the usual supergravity scenarios (where 
there is no problem), the extension of the theory by adding a weak-scale singlet is 
more justifiable, because it circumvents the unwanted relation B^/ /i ~ F/M, it does not add 
many new free parameters, and it cannot destabilize the hierarchy, at least for sufficiently 
low mediation scale. 

In this paper, we have studied a model of gauge- mediation with the addition of a singlet, 
including a messenger-singlet coupling, that was first proposed in ref. [23] (we denote this 
model as N-GMSB). We have computed the induced soft terms and analyzed the phenomeno- 
logical consequences. The singlet interactions are described by the three couplings A, k and 
^jj. However, since two of them can be traded for /x and B^, the theory contains only one 
extra parameter with respect to the ordinary gauge-mediated supersymmetric SM without 
singlet, and therefore maintains its high predictive power. Electroweak breaking requires a 
mild fine tuning, endemic to supersymmetric models after LEP2. However, the necessary 
"critical" condition can be achieved for a large range of values of the coupling constants A, 
k and (at the price of tuning one of the three parameters). 

The lightest Higgs mass gives the most stringent constraint to N-GMSB, leading to a 
heavy supersymmetric mass spectrum and a large messenger scale M, and identifying three 
special regions in the space of couplings A, k and ^u. Region I is characterized by large 
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tan (3 and light singlet fields. The singlet pseudoscalar is particularly light and the singhno 
can be the NLSP, leading to a potentially characteristic signal of the supersymmetric decay 
chains, with Higgs and missing energy in the final states. Region II has large tan /3, heavy 
singlet fields and a low-energy mass spectrum that is very similar to the one of usual gauge 
mediation. Region III has tan (3 close to one, and the tree-level quartic Higgs coupling is 
mostly generated by A. Anyway, even in these regions the lightest Higgs mass is not larger 
than the maximal value attainable in the usual GMSB. We find that the Higgs can be 
substantially heavier than 120 GeV only if we abandon perturbativity of the couplings up 
to the GUT scale. 
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Appendix A 



In this appendix we provide the RG equations for the gauge and superpotential couphngs of 
N-GMSB, vahd above the messenger scale M. The RG equations for the couphngs and the 
soft supersymmetry-breaking terms of the N-MSSM, vahd below the scale M, can be found 
e.g. in ref. [22] (note however that our definition of A\ differs by a sign from that of ref. [22]). 

Defining [3\^ as in eq. (ITUl) . the RG equations for the couplings Aj are: 

^ = 1^(^ + 1)' (48) 
^s = 1^(^-3), (49) 
Ph. = ^ (e/^? + hi + X' -3g' - fg- - f gA , (50) 



= ^jQhl + hl + hl + X'-3g^-lg'^~fg^,), (51) 



Ph. = ^A^K + 3hl + X'-3g'-3g"] , (52) 

/3, = :^{AX^ + 2e + 3h^t+3hl + hl + 2eD + 3eT-^9'~9"] , (53) 

Pk = |^(6A^ + 6A;^ + 6^^ + 9e^) , (54) 

P,, = ^(2A2 + 2P + 4e!, + 3e^-3/-^'^) , (55) 



2y + 2k' + 2eo + ^eT-^9i-7;g"] , (56) 



where n is the number of messenger pairs. 
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Appendix B 



In this appendix we derive the expressions for the soft supersymmetry breaking terms in the 
scalar sector using the wave-function renormahzation method proposed in refs. 



The soft terms in eqs. f|TT|) - f|T2|) are given in terms of the following derivatives of the 
wave-function renormahzation Z with respect to the messenger mass M, evaluated at the 
renormahzation scale Q (for simplicity we will drop the field index a in this appendix) 



glnZ(M,g) ^„ ^dHxvZ{M,Q) 



2 (91nM 



4 (9(lnM)2 



(57) 



By integrating eq. ([9]) between an arbitrary high-energy scale A and the renormahzation 
scale Q (with Q < M), we obtain 



In 



Z{Q) 
Z{A) 



rlnM flnQ 

/ dt 7^+) + / 

Jin A JlnM 



dt 7^-) 



(58) 



where 7*^^-* are the anomalous dimensions above and below the messenger scale M, respec- 
tively. 

Taking the first derivative of eq. (!58l) . we obtain the expression of Z'{Q = M) shown in 
eq. (|T3l) . Taking the second derivative, we find 



Z"\ 



1 



a(A7)aAf(M) 97(") 9Af(g) 



d\J d\nM d\f d\nM 



=M 



(59) 



With the help of eq. (fTOl) . we obtain 

9A^(M) _ (+), 
d\nM I 



=*^' d\nM 



(60) 



Replacing eq. fl60|) into eq. fl59|) . we obtain the expression of Z"{Q = M) shown in eq. ([T3 
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